Yes, this looks like the right direction. We only need the vertical sum of AOD though, not every layer. And I think if you could output AOD in this test: https://github.com/CliMA/RRTMGP.jl/blob/main/test/all_sky_with_aerosols_utils.jl, it would be a good step towards having it in ClimaAtmos.

@sriharshakandala, can you please take a look at this?

pseudocode

    band_wv_550nm = 1 / band_550nm
    bandidx = 0
    for i in 1:nband
        if bnd_lims_wn[1, i] <= band_wd_550nm <= bnd_lims_wn[2, i]
            bandidx = i
            break
         end
    end

I'm fine with incrementing the minor version in this PR, too. cc @szy21

I'm fine with incrementing the minor version in this PR, too. cc @szy21

Oh, right, I forgot to do it.